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Introduction 


The main objective of this research is to obtain a calibration 
of a low concentration of charged defects in some III-V semiconducting 

compounds. The experimental technique being used is nuclear magnetic 

/ 

resonance (NMR), and the objective is to be attained through a thorough 
examination of the properties of the NMR lines of the nuclei in the III-V 
compounds . 

Understanding the properties of the NMR lines, and how they are 
influenced by the presence of charged defects, requires both theoretical 
and experimental exploration. This report is thus divided into three sec- 
tions: I. Theory, II. Experiment , and III. Conclusions. 



I . THEORY 


The sensitivity to defects achieved in this work will depend on a 

detailed knowledge of the magnetic resonance line shape. The best calculation 

of the magnetic resonance line shape in solids was done by Lowe and Norberg 

( 2 ) 

However Abragam' has proposed a two parameter shape function that reproduces 

the observed shapes as well as the Lowe-Norberg theory and is much simpler in 

( 3 ) 

form. More recently Betsuyaku has extended the Lowe-Norberg calculation to 
include higher order terms in their expansion. This extension removes the agree 
ment between the theory and experiment. Thus this fundamental problem remains 
in unsatisfactory shape. 

A new approach to the problem is needed. It is suggested that a 
combination of "single site approximations" and the "coherent potential approxi- 
mation"^ may prove useful. 

In order to apply the methods suggested above, a knowledge of the 
normal modes, or at least a set of almost normal modes, of the spin system is 
needed. A preliminary and still incomplete study of these modes is described 
below. 

The secular part of the dipole-dipole Hamiltonian coupling a 
system of identical spins is 
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y is the magnetogyric ratio for the spins, 
ft is Planck's constant, 

T is the lattice vector designating a nuclear position, 

_ T-f 

iWJi 

/-} o is the external d.c. magnetic field, 

is the spin operator of the nucleus, and 
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The objective is to find a normal modes transformation that diagonalizes 
this Hamiltonian. 

Try the transformation 




it-X 


(3) 


where N is the number of spins in the system, including periodic boundary 
conditions, it is a reciprocal lattice vector which ranges over the first 
Brillouin zone, and ^ is a polarization index which takes on three possi- 
ble values. As the calculation proceeds it will become obvious how fa 


and 




are to be found. 


Substitute Eq. (3) into Eq. (l). Then one finds 
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where 


Mr) h n i™'-* -- M-r) 
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Choose the €+>/ to be the eigenvalues of the diadic <A0, i-e., 

= * w V' ( s' %'$ ■ 
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Since ($ / ) rr t/L (~&% we know that 


and 


cA is a diagonal diadic, hence its eigenvalues 
and eigenvectors are trivial to find. This general formalism is presented in 


k-Vp' = ~ 


case it ever becomes useful to use this method on more complex spin Hamiltoni- 


ans, e.g. the full dipole-dipole interaction. The eigenvectors of 


JUt) 


are 
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and the corresponding eigenvalues are 
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with 005(9/ */y o . Note that 60 0 ^ r 0 for a crystal with 


center of symmetry. However we shall see that UJg ^ is almost discontinuous 

near k. = O and, for k ~ ^ is not small. An approximate 

expression for 60 is developed in Appendix I. 
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Nov return to the diagonalization of • Inserting Eqs. (5) and 
(6) into Eq. (4) converts it to the form 


" A/ 4^ 

The sum over X 1S 
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2_, e - Og -£/ , a Kronecker delta of K. 


and tc. Then, since ^ ^ ^ , the becomes 
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Since 1-* is a Hermitean operator, Eq. (3) gives us 
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where the last equality is reached by substituting -tc for tc everywhere 

^ r\ 

in the term on the left and using the fact that e t?(3 " e -tf ‘ ThS 

equality of the last two terms in Eq. (10) implies that 


f-Sf 1 > 


an 


since the expressions must hold for arbitrary X • Eq. (ll) can be used 
to write Eq. (9) in the form 
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In order to complete this discussion, we now examine the commuta- 
tion relations among the f* ^ and °P erators * First the inverse 

transform of Eq. (3) yields 


fa® = h £** ' ? 


-if.je 


(13) 


Then the commutator between 


is 


but 
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Examine the special cases (3; = x >3> ? : 


c W>fa'^ s - 0 ' 

, <-A ,-tA -tCfe+fc^’X t 
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Similarly, one finds 


0 ) 

i 

tytt* 4*3^ ft > 

uv4y =0) 

[■ ijtx » ffe'J = w , 

plus cyclic permutations of Eqs . (20) and (22). 

Using Eq. (19), Eq. (12) can be rewritten in a familiar form: 


• (23> 

Although this Hamiltonian is written as a sum of one "particle" 
operators, the operators corresponding to different polarizations do 
not commute. Hence the Hamiltonian is not yet diagonal. Operators for 
different but the same polarization do commute. There are a few other 
commutation relations and theorems which can be derived, but so far we have 
been unable to find another transformation that completes the diagonalization 
of the Hamiltonian. 

For example, the net magnetization can be expressed in terms 

of the as 

= £ 4 1 ,? = • 


( 24 ) 
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Other commutation relations that prove to be useful are (in these expressions 
the arrow is left off the vectors , Jc ): 



From Eqs . (25), (26 ), and the fact that 
is easy to show that 






it 


< 31 > 

or equivalently from Eq. (24) 
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Eq. (32) is not too surprising since that is the condition used to pick the 
secular part of the full dipole-dipole Hamiltonian. 

if X can be completely diagonalized then the Lowe-Norberg pro- 
blem can be solved exactly. However, since this Hamiltonian is a generaliza- 
tion of the Ising model which has yet to be solved, it seems unlikely that we 
shall succeed. It does seem worth a bit more effort. The fact that the 
commutation relations, Eqs. (l5)-(22), involve a factor may prove to 

be useful in some expansions. However care must be used since some 
may be of order >/aT . For example, for a system of spin 1/2 particles in 
the ferromagnetic state with all spins up lt> 




A 

— 2_ n/kJ kl_v « 0 


X 
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the k- - 0 matrix element is large while all the others are zero. For states 
with nearly random spin orientations, the matrix elements of the >S 

are probably mostly of order , but this hasn't been proven. 

A good deal remains to be done on this subject so this work will 

( 4 ) 

continue. We shall also try to extend these ideas to help analyze the Waugh 
type experiments. 



II . EXPERIMENT 


E. H. Rhoderick has investigated the nuclear magnetic resonance of 
doped and undoped indium antimonide and gallium arsenide. ^ He found that the 
intensity of the resonance absorption derivative decreased with increasing con- 
centrations of impurity. In fact, his data show a sharp decrease in the inten- 

3. 3 3 

sity of the resonance as the impurity concentration increases up to 10 cm 

1 8 3 

Concentrations larger than 10 cm serve to further reduce the intensity, but 
at a much slower rate. 

The explanation for this behavior is as follows. In the pure com- 
pound the resonance line is broadened significantly mostly by the magnetic 
dipole-dipole interaction. The resonance line in a perfect crystal is not 
broadened by the quadrupole interaction because the. lattice has cubic symmetry 
and each atomic site has zero electric field gradient. An impurity destroys the 
cubic symmetry in its neighborhood, and broadens the resonance line through the 
quadrupole interaction. Since the quadrupole moments of each of these nuclei 
are large, the satellite transitions (3/2 -*►1/2, etc.) are split far from the 
resonance line and do not contribute to the detected resonance. The initial 
sharp decrease in intensity is thus due to the loss of the satellites from the 

resonance line. The gradual decrease in intensity observed by increasing impurity 

10 3 

concentrations greater than 10 cm J is due to a broadening of the central com- 
ponent (the 1/2 -*>1/2 transition) by the quadrupole interaction. This broadening 
is always present, but its effect can be distinguished only when the central 
component alone is detected. 

Among the other conclusions that Rhoderick draws from his data is that 

the antishielding factor necessary to explain the sharp decrease in intensity 

3 

is very large, on the order of 10 . We have used Rhoderick' s data to make an 
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- 10 


independent calculation of the strength of the quadrupole interaction, and we 
find that it is 2000 times stronger than predicted by a simple point charge 
defect model. 

Our calculation consists of assuming that the impurity atom will assume 
a net charge at its lattice site of r e , where ^ is a dimensionless parameter 
describing the strength of the interaction. This produces an electric field 
gradient at the sites of nearby nuclei. In particular, a nucleus a distance r 
away from the impurity will have satellite lines shifted by an amount 



where £ is the dielectric constant, cQ is the electric quadrupole moment of 

the nucleus, and 6 is the angle between the vector from the impurity to the 

nucleus and the static magnetic field direction. If f is small, the satellite 

shift will serve only to broaden the resonance line; but a large f g corresponds 

to loss of the satellite components of the resonance. In this calculation of 

the strength of the quadrupole interaction we choose a frequency f such that if 

f f for a given nucleus, then we may still see its satellite contributions; 
s 

but if f > f, it will be split far enough that we cannot detect it. 

s *" 

In order to calculate the total number of nuclei that contribute only 

the (-1/2 -*►1/2) component to the resonance line we need to calculate the volume 

69 

about each impurity in which f > f for every Ga atom. This volume is found 

s 

to be , 

,kflV3/Xe_Q \/££±3\ 

V = TfW hf s ]( -re ' • 

Then the number of nuclei within the volume is N = n Q V, where n Q is the atomic 
69 

density of Ga . Since, as we shall see, this interaction is very strong, then 



XI 


the volume V is large and we are justified in using a continuum model for the 
lattice . 

The intensity of the resonance line is proportional to the density of 
Ga^ atoms, n , within the specimen; but the intensity of the resonance when 
all the satellites are shifted beyond detection must be proportional to .k n Q . 

If the interpretation of the intensity behavior is correct, then the linear de- 
crease of intensity with increasing impurity concentration should extrapolate 
back at zero concentration to the intensity of an unbroadened central component, 
which is ,1+n . Actually, Rhoderick's data extrapolate back to .6 of initial in- 
tensity, indicating that in his "pure" crystal not all the satellites were present. 
Let n' represent the atomic density which at zero impurity concentration contri- 
butes only the central component . Then with zero impurity concentration the in- 
tensity may be written 

where K is the normalization constant. Rhoderick normalizes this intensity to 
unity. The extrapolation to zero concentration yields 

. u = .(o K (n b - kn') = .4-K^o } 

or 

, y 3 

r) b -~.kr> - 5^0 } K=2 .h 0 

Rhoderick's "pure" sample contained only k/9 of the satellites. We can now write 


the intensity of the resonance at a given impurity concentration, C, 


and by differentiation. 


I = (n 6 -.fcn' — .tn 0 V c), 




- \~ts- ' 1 

We use the value of from Rhoderick's data and the value of 3 KHz 

go 

for f , which is half the total width of the Ga resonance line, and obtain the 

S * 
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result ^ = 2000. Since the impurities used by Rhoderick were one electron 
donors or acceptors, we are reasonably confident that the effective charge 
at the impurity site is nearly + e. Thus the effective strength of the 
quadrupole interaction is a factor of 2000 larger than predicted by this 
calculation for a point charge in a continuous dielectric medium. 

We have made preliminary measurements on a GaAs sample that show 
qualitatively the same results as Rhoderick. In our case, however, we 
have induced vacancies in the sample by thermal damage - heating the sample 
to about 700°C, then rapidly quenching it to about room temperature. These 
vacancies then act as charge centers and create electric field gradients 
like impurities. We have observed both a decrease in NMR signal intensity 
and a change in NMR linewidth. 

The decrease in signal is shown in Fig. 1 for the one sample ob- 
served so far. There is an orientation dependence of the signal intensity, 
which appears to be preserved even after the damage. Whether or not this 
is a real effect cannot be determined yet, since these numbers were obtained 
by measuring the peak-to-peak height of the derivative of the absorption 
curve , X". We use this number here because Rhoderick did, and we were 

interested in comparison with his data. However, this measurement does 
not accurately describe the changes in %" as impurities are introduced 
into the lattice, for if changes shape as well as intensity the single 


measurement of 


ax" 

an 


J mi, 

4m V 


might well lead to an incorrect interpretation of 

'max 


the n umb er of nuclei contributing to the resonance. 


We feel that it is incorrect to attribute the enhanced quadrupole 


effects caused by lattice defects - either impurities or vacancies - to 
the ordinary antishielding action of the electrons. In fact there is some 



- 13 - 


indication that another process is involved, which is not totally under- 
stood, and which may not have been seen by Rhoderick. In Rhoderick's 
"pure" samples he found no line width change as a function of orientation 
in the magnetic field, and he explains this isotropic broadening using the 
pseudo-exchange interaction, in which the nuclei interact with each other 
through electrons . 

( 6 ) 

R. K. Sundfors has used HMR and nuclear acoustic resonance 

techniques to isolate broadening mechanisms in III-V compounds. Unlike 

69 

Rhoderick he finds a dependence of the Ga second moment upon crystal 

orientation in the magnetic field. He attributes this to quadrupole 

broadening caused by a small random electric field gradient present in 

the lattice. We find the same angular dependence in an undoped specimen 
5 

with 10 fewer charge carriers than Sundfors'. We feel that fewer charge 
carriers implies fewer substitutional impurities, which should cause much 
less broadening due to random electric field gradients. When damage is 
caused by thermal shock we find that the linewidth decreases in addition 
to the loss of intensity. These line width data are shown in Fig. 2. This 
behavior indicates that spin-spin relaxation is partially inhibited when 
the quadrupolar satellite shift occurs. Thus the sharp decrease in "in- 
tensity" (as defined by Rhoderick) may well not be entirely due to anti- 
shielding. 

Some of these data were presented in an oral report at an American 

( 7 ) 

Physical Society meeting, and the abstract for that report is included 
here as Appendix II . 

The best way to monitor intensity is to integrate over the ab- 
sorption curve , but there is a difficulty in that a low signal to 



noise ratio in the extreme wings of the resonance causes very large experi- 
mental error. In order to make such measurements accurately, one must use 
signal averaging techniques. Much of the work of this past year has been 
devoted to purchasing and building the components of a system to do this. 
The system consists of an analog-to-digital converter, a small computer, 
and the interfacing required to make the various parts of the experiment 
work together. The signal averaging system is nearly complete, and is 
expected to be operating about June 1, 1972. 

A pulsed NMR spectrometer has been completed and is now operating 

This apparatus will assist in careful line shape and intensity measurements 

(I*) 

and as well allow the application of the Waugh pulse techniques to the 


problem. 



III. CONCLUSIONS 


Substantial progress has been made in preparation for the experi- 
mental and theoretical research on the main objective of this contract. 

The preliminary investigations have already found new and interesting re- 
sults, some of which are being published. The sensitivity to thermally 
induced vacancies that we have found supports our hypothesis that we will 
be able to calibrate low concentrations of charged defects. A quantitative 
evaluation of the ultimate sensitivity of this technique must await more 
detailed studies of the NMR line shape and intensity . This phase of the 
experimental program will soon be aided considerably by the recently acquired 
signal averaging system. 

It was hoped that the signal averaging system would be completed 
by the time of this report. The major component of that system is the small 
digital computer, for which an order was first place in late July, 1971. 
However, the order was held up by the Division of Automatic Data Processing 
under the Virginia Governor's office until November 18, 1971. This delay 
was purportedly to coordinate our purchase with state acquisition of other 
computers, although this federally- funded research project was not intended 
to be coordinated with other state computers or research projects. The 
outcome of this is a delay of approximately four months in beginning our 
full effort in taking data. 

Some of the theoretical research done this past summer was carried 
out by A. Sher at the Aspen Summer Institute for Physics in Aspen, Colorado. 
This was done with the knowledge and consent of both the principal investi- 
gator and the technical officer of this contract. 
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Appendix I 


To gain some idea of the behavior of » convert the sum over 

to an integral from Ji - GL to fL - oo . If the upper limit on /- 
is taken to be the size of the sample f\l 3 , then is zero; if jL~°° 

is used, then (D ox is finite. Since this is a single point in k -space, 
it doesn't contribute much to integrals over 1c which will arise later in 
the theory; so practically no error is introduced. Also assume that the 
radial distribution function is uniform and equals '/a. 3 , i .e. , this is a 
first approximation to a simple cubic lattice. 

Denote a spherical harmonic by y:oi> , where L. is the angular 
momentum , M the 2 component of the angular momentum, and - © j ^ . 

The second spherical harmonic Y,°cn) is 


T CrO = (fet ( 3 - 1} • 

Z 


(A-l) 


The function £. 


I £•! 


can be expanded in spherical harmonics as 


e l * J =4T r(&)t i L Y;W>^<w) (A. 

L.fA 


-2 ) 


where o; (x) is a Bessel function of order V with argument X 

Using Eqs . (A-l) and (A-2), and converting the sum in Eq. (7) to an integral 
S- J ^3 dT = J ■£/(£ dX dSL , the expression for 
becomes 
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But the integral is 


JjaY; - V, S ( 


and when this is inserted into Eq. (A-3)» it reduces to 


^ (-fort ^ _ 

a? 

The integral in Eq. ( A-b ) over JL equals k^'J} CkdS) /(I c(l) 1 and 
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T /y 1 ! — ( ^ - C ^5- X 1 so Eq. (A-4) can he simplified to 
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The frequency dispersion has the form 




Appendix II 


Abstract Submitted 
for the Atlantic City Meeting of the 
American Physical Society 


March 27-30 , 1972 

Published in: Bull. Am. Phys . Soc. 17., 332 (1972). 


Physical Review 
Analytic Subject Index 
Number 21 . 4 


Bulletin Subject Heading 
in which Paper should be placed 
Nuclear Magnetic Resonance 


Anti shielding Factor From Gallium Arsenide NMR . * 

JOES? , College of William and 


R. HESTER, A. SHER , J. P 
Mary — Using the data, of 
variation 


H. Rhoderick- 1 - on the intensity 
-che Ga^9 MR absorption derivative line 
with a change in the Te concentration of Te-doped GaAs , 


we 


calculate an effective antishielding factor on tne 

lieve that the Ga°9 


order of 1000. This leads us oO oe 
resonance can be sensitive enough to detect defect con- 
centrations in undoped GaAs possibly as small as 10 ' - 
1015 cm -3 . Preliminary data taken by cv NMR have shown 
qualitatively the same results as Rhoderick has presented. 
The creation of defects by thermal damagp, in the same 
manner as H. R. Potts and G. L. Pearson 2 , lowers the in- 
tensity of the resonance line, ana. m audition, causes a 
line narrowing. This line narrowing may be the result oi 
a partial quenching of one or more of the spin-spin 
interactions. 


^'Supported in part by NASA. 

1r. H. Rhoderick, J-. Phys. Chem. Sol. 8_, 498 *1958). 

2 H. R. Potts and G. L. Pearson, J. Appl. Phys. 47., 2098 

( 1966 ). 
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